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Abstract 

The  temperature  profile  inside  a  superconducting  target  in  laser 
ablation  is  calculated  for  laser  pulses  of  various  shapes.  The  calculation 
is  based  on  the  equation  of  heat  conduction.  All  parameters  characterizing 
the  target  material  are  assumed  to  be  temperature  dependent  and  are 
determined  empirically  by  extrapolating  experimental  data  to  the  melting 
point.  The  receding  velocity  of  the  vapor -solid  interface  is  determined  by 
the  dynamical  balance  of  energy.  Our  calculation  shows  that,  in  general, 
there  exist  subsurface  overheating  spots  at  different  instants  as  long  as 
the  laser  pulse  intensity  is  strong  enough.  The  dependence  of  their 
occurrence  on  the  pulse  shape  is  analyzed,  and  conditions  to  avoid  them 
without  jeopardizing  the  deposition  process  are  discussed.  Ac9,-» ■> tou  r 
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I. 


Introduction 


The  discovery  of  high-Tc  superconductivity  has  prompted  numerous 
efforts  in  the  preparation  of  high-quality  superconducting  thin  films. 


Among  many  successful  techniques  to  date,  such  as  sputtering,  thermal  and 
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electron  beam  evaporation,  molecular  deposition  and  sol-gel,  the  pulsed 

5 IL 

laser  deposition  (PLD)  ^ its  unique  advantages.  The  PLD  technique  is 
relatively  inexpensive,  and  provides  high  deposition  rates  and  excellent 
control  of  ^stoichiometry  of  the  deposited  film.  It  even  eliminates  the 
post -annealing  procedure.  Furthermore,  the  low  chamber  pressure  during  the 
deposition  process  reduces  the  undesirable  impurities. 

Much  progress  has  been  made  recently.  For  example,  Y-Ba-Cu-0/Yx-Pr^_x- 


Ba-Cu-0  superlattices  have  been  successfully  grown  by  the  PLD  technique. ** 
Techniques  are  also  greatly  advanced  in  the  fabrication  of  controllable 

Josephson  junctions^  and  in  the  growth  of  thin  films  of  YBaCuO  (YBCO)  onto 

O 

the  new  magnetic  substrate  YbFeO-j .  The  chamber  temperature  during  the  PLD 


process,  however,  has  to  be  greatly  reduced  if  the  technique  is  applied  in 
the  integration  of  superconductors  and  semiconductors  for  device  purposes. 
This  is  because  most  semiconductors  can  not  tolerate  temperatures  higher 
than  450aC.  The  introduction  of  an  oxygen  jet  into  the  laser-ablated  plasma 
has  made  it  possible  to  process  superconducting  thin  films  at  650°C. 

Further  reductions  in  temperature  have  been  reported  by  the  incorporation  of 
9 

an  oxygen  plasma.  A  new  technique  involving  plasma-assisted  laser 
deposition  has  been  developed  to  fabricate  superconducting  YBCO  films 

without  post-annealing.^ 
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The  quality  of  superconducting  films  made  in  the  PLD  process  depends  on 
a  number  of  factors  such  as  the  energy  density,  wavelength,  duration  and 
shape  of  the  laser  pulse,  the  incident  angle,  the  optical  and  thermal 
properties  of  the  pellet,  and  the  distance  between  the  pellet  and  the 
substrate.  As  a  matter  of  fact,  the  PLD  process  consists  of  three  stages: 

1)  laser  evaporation  (LE)  of  the  superconducting  target,  2)  transport  of  the 
evaporated  material  from  the  target  to  the  substrate  and  3)  deposition  on 
the  substrate.  The  evaporated  material  may  be  atoms,  ions,  molecules  and 
clusters,  and  the  presence  of  large  molecules  and  clusters  eventually 
influences  the  film  quality.  The  transport  process  is  most  complicated  to 
control  experimentally  and  most  difficult  to  treat  theoretically.  Particles 
of  various  sizes  and  charges  move  in  an  oxygen- filled  chamber  exposed  to 
laser  pulses,  and  hence  collisions  of  all  types  are  involved  in  this  stage. 
The  final  stage  involves  the  deposition  of  all  particles  including  products 
of  chemical  reactions  during  the  transport. 

Attempts  have  been  made  to  understand  the  nature  of  these  processes. 

The  angular  distribution  of  the  composition  and  thickness  of  the  YBCO  film 

has  been  measured. ^  It  appears  that  there  exist  two  distinct  components  in 
the  deposited  material.  One  is  the  nonstoichiometric  cosff  component,  where 
9  is  the  angle  with  respect  to  the  normal  to  the  substrate  surface,  and  the 
other  is  the  stoichiometric  forward- directed  component.  An  anisotropic 
expansion  of  the  high-pressure  laser-evaporated  plasma  has  been  suggested  to 

12 

explain  such  spatial  variations  in  the  deposition.  During  the  LE  process, 

13 

a  partially- ionized  plasma  has  been  observed  above  the  target  surface,  and 
high-energy  atomic  beams^  and  cluster  emmision^  have  also  been  found  in 
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experiments.  Besides,  it  is  known  experimentally  that  the  pulse  duration 

plays  an  important  role  in  the  final  film  quality. ^ 

On  the  other  hand,  much  less  theoretical  research  has  been  conducted  to 
understand  the  complex  mechanism  governing  the  PLD.  A  model  was  developed 

with  predictions  in  qualitative  agreement  with  experiments.^  It  assumes 
that  a  high- temperature  plasma  is  created  out  of  the  evaporated  material 
under  the  irradiation  by  the  laser  pulse.  An  adiabatic  expansion  of  the 
plasma  then  results  in  the  deposition.  The  treatment  of  the  first  stage  of 
LE,  however,  is  oversimplified  in  this  model,  so  that  it  fails  to  provide 
the  temperature  profile  of  the  target.  The  change  of  temperature  with  time 
and  depth  in  the  target  is  necessary  to  investigate  the  possibility  of 
18 

explosion,  since  liquid  drops  as  well  as  solid  particles  have  been 
19 

observed. 

20 

A  recent  calculation  of  the  temperature  profile  in  the  high-Tc 

superconductor  target  irradiated  by  laser  pulses  of  triangular  shape  shows 
no  subsurface  thermal  peak  which  is  crucial  to  the  subsurface  explosion. 

The  target  surface  is,  however,  simply  specified  in  Ref.  20  by  the  melting 
temperature,  and  all  parameters  are  assumed  to  be  temperature  independent. 

21 

A  more  careful  calculation  shows  the  existence  of  subsurface  overheating 
spots.  On  the  other  hand,  subsurface  thermal  peaks  do  exist  and  have  been 

22  23 

observed  in  other  materials.  ’  There  is  certainly  no  evidence  showing 
why  the  subsurface  explosion  can  not  occur  in  high-Tc  superconductors. 

We  consider  in  this  paper  a  YBCO  target  exposed  to  laser  pulses  of 
different  shapes  and  intensities,  as  shown  in  Fig.  1,  and  study  its 
evaporation.  We  attempt  to  calculate  the  temperature  profile  inside  the 
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material  and  its  variation  with  the  pulse  shape  and  intensity.  As  the 
target  surface  is  receding  during  the  evaporation  process,  we  solve  the 
equations  in  the  frame  that  moves  with  the  surface.  Temperature - dependent 
parameters  are  determined  by  extrapolation  to  high  temperatures  from  the 
empirical  formula  constructed  on  the  basis  of  experimental  data. 

II.  Method  of  calculation 

When  the  superconducting  target  material  is  irradiated  by  an  intense 
laser  pulse,  its  surface  layers  are  heated  continuously  until  the  outer 
surface  reaches  the  vaporization  temperature  T  .  Then  the  material  is 

removed  from  the  target  surface  in  the  form  of  electrons  and  positive  ions, 
molecular  clusters  or  solid  particles.  Due  to  further  absorption  of  energy 
from  the  laser,  the  energy  absorbed  by  the  target  forms  a  distributed  heat 
source  around  its  surface  region,  and  this  region  is  mobile  because  of  the 
laser  evaporation.  Our  problem  is  to  study  the  motion  of  this  vapor-solid 
18  20  23 

interface  ’  ’  and  the  temperature  distribution  into  the  target  material. 

During  the  laser  irradiation,  the  temperature  T(x,t)  in  the  material  is 
a  function  of  both  the  time  t  and  position  x  >  X,  where  X  is  the  depth  from 
which  material  has  been  removed  by  the  laser.  The  lateral  rate  of  heat 
diffusion  is  assumed  to  be  small  in  our  treatment.  This  function  is 
determined  by  the  heat  diffusion  equation 

2 

pm  Cp(T)  J ~  T(x,  t)  -  a(T)b(T)I<t)e'b(T)(3C'X)  +  *<T)  ^  T(x,t),  (1) 

3xZ 

where  we  have  assumed  the  most  general  temperature -dependent  parameters, 
namely  the  density  p,  specific  heat  Cp,  thermal  conductivity  k  and 
absorption  coefficient  b.  I  is  the  time -dependent  laser  intensity,  and  the 
function  a(T)  is  introduced  to  account  for  the  reflection  loss  of  the 
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surface  and  the  screening  effect  of  the  vaporized  material.  This  is 
necessary  because  the  vaporized  material  from  the  target  absorbs  part  of  the 
laser  energy  and  hence  reduces  the  energy  the  target  can  absorb.  The  energy 
balance  at  the  interface  requires  that  the  energy  given  to  the  vaporized 
material  is  equal  to  the  heat  energy  flow  to  the  surface,  that  is, 

XWV  |f  -  *<T>  fc  T(x,  t)j  X_x  .  (2) 

where  Lv  is  the  vaporization  heat. 

The  temperature -dependent  parameters  plus  time -dependence  of  the  laser 
intensity  make  it  extremely  difficult,  if  not  impossible,  to  solve  Eqs.  (1) 
and  (2)  analytically.  Before  we  proceed  to  solve  these  equations 
numerically,  we  first  transform  them  into  dimensionless  forms  by  introducing 
the  following  variables:  6  -  T/Tv,  r  -  t/t q,  S  -  (x-X)/iQ,  L  -  pLv,  C  -  pCp 

and  B  -  Ibtg/CT^.  We  have  chosen  in  these  definitions  tQ  -  30  ns,  a  typical 
laser  pulse  duration  and  fg  -  3  x  10'^cm.  The  laser  intensity  is  measured 


8  2 

in  the  unit  of  1.27  x  10  W/cm  .  The  parameter  B  is  introduced  just  for 


convenience.  Thus,  Eqs.  (1)  and  (2)  become 
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where  u  stands  for  the  speed  of  the  moving  interface.  Equation  (4) 
indicates  the  fact  that  the  interface  does  not  move  until  the  surface 
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material  at  high  temperatures,  however,  have  never  been  measured,  at  least 
to  the  best  of  our  knowledge.  Therefore  we  have  to  determine  the 
temperature  dependence  of  the  parameters  empirically  from  the  available  data 
before  we  can  carry  out  the  numerical  work. 

For  the  1:2:3  material,  we  construct  the  approximate  formulas  from 

24-28 

large  amount  of  experimental  data  by  the  extrapolation  method. 

According  to  the  experimental  curves  given  in  Ref.  24,  we  find  by 
extrapolation  that  the  thermal  conductivity  can  approximately  represented  by 
the  formula 

k  -  0.042  +  A  AT  (5) 


in  the  unit  of  W/K-cm,  where  AT  —  T  -  300K  and  A  —  1.2  x  10  ^W/K^cm.  The 
density  is  a  linear  function  of  T  according  to  Ref.  25.  The  determination 
of  the  specific  heat  is  more  complicated.  From  the  thermal  expansion 

27 

coefficients  in  Ref.  26  and  the  measured  data  of  the  density,  we  can 
deduce  the  relation  between  p  and  T.  The  temperature  dependence  of  Cp  can 
be  deduced  from  Ref .  28 ,  and  finally  we  find 

C  -  1.82  +  2.6  x  10‘8  AT  (6) 


in  the  unit  of  J/g-K.  The  absorption  coefficient  is  a  result  of 
extrapolation  of  the  data  presented  in  Ref.  27,  namely, 

b  -  2.27  x  105(1  •  9  x  106  AT)  (7) 
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in  the  unit  of  cm"^.  For  the  screening  factor  we  take  a  -  0.9  for  6  <  1  or 
before  the  vaporization  starts,  and  a  -  0.3  when  the  vaporization  begins  ($ 

29 

-  1).  This  choice  agrees  with  the  experiment.  It  is  noted  that  the 

20 

boiling  temperature  is  still  not  known.  A  rough  estimate  is  twice  the 

-an 

melting  temperature,  which  is  about  1300K.  We  take  -  3000K  in  our 

31 

computation.  We  have  in  fact  found  that  a  change  of  Tv  within  the  range 
of  2000  -  3000K  does  not  change  the  result  qualitatively. 

III.  Results  and  discussion 

The  temperature  change  with  the  depth  into  the  material  is  calculated 
for  various  intensities  of  incident  laser  pulses.  For  all  pulse  shapes  we 
have  considered,  it  is  possible  to  find  subsurface  temperatures  higher  than 
the  surface  temperature  for  sufficiently  high  intensity  of  the  pulsed  laser. 
This  means  that  internal  explosions  are  possible  and  hence  provide^  a 
reasonable  explanation  of  the  observed  liquid  droplets  and  small  solid 
19 

particles.  The  subsurface  overheating  phenomenon  can  be  understood  as  the 
penetration  of  the  laser  field  into  the  material  and  the  creation  of  heat 
sources  therein.  On  the  other  hand,  the  rate  of  heat  dissipation  on  the 
surface  is  much  greater  than  that  in  the  interior  because  of  the  surface 
evaporation. 

Temperature  profiles  calculated  for  different  cases  are  plotted  in 
Figs.  2-5,  and  the  corresponding  interface  receding  velocity  in  each  case  is 
listed  as  a  function  of  time  in  Table  1.  For  a  rectangular  pulse,  the  laser 
intensity  remains  constant.  The  interface  receding  acceleration  increases 
at  first  and  then  decreases,  and  the  subsurface  peak  temperature  also  shows 
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saturation.  This  is  because  the  rate  of  thermal  diffusion  is  proportional 
to  the  temperature  gradiant.  The  situation  is  depicted  in  Fig.  2(a)  and 
Table  1(a).  Figure  2(b)  and  Table  2(b)  show  the  results  for  a  triangular 
pulse.  As  the  laser  intensity  increases  with  time  until  its  peak,  the 
interface  receding  speed  also  increases.  As  a  result,  the  subsurface 
thermal  peak  moves  into  the  material  continuously  with  increasing  peak 
temperature.  In  other  words,  the  saturation  phenomena  are  not  apparent 
because  the  laser  intensity  keeps  increasing. 

To  see  how  the  time  spread  of  laser  pulses  of  the  same  peak  intensity 
affect  the  temperature  distribution,  we  take  the  triangular  pulse  as  an 
example.  We  consider  an  ideal  case  in  which  the  pulse  heats  the  interface 
up  to  its  vaporization  temperature  T  as  the  intensity  increases  linearly 

with  time  and  reaches  its  peak.  Now  we  let  the  pulse  intensity  decrease  at 
different  rates  by  assuming  d  -  1  and  0.17  in  Fig.  1(b).  The  results  are 
plotted  in  Fig.  3.  It  is  observed  that  in  Fig.  3(a)  and  Table  1(c)  the 
subsurface  thermal  peak  still  exists.  The  peak  temperature  (0  -  1.08)  is 
considerably  lower  than  is  Fig.  2(a)  ( 8  -  1.7)  because  the  laser  peak 
intensity  has  been  reduced  from  unity  to  0.28.  The  interesting  point  is 
that  even  if  the  laser  intensity  decreases  linearly  with  time,  the  peak 
temperature  is  still  increasing  with  time  during  0  <  r  <  0.24  and  the  peak 

position  keeps  moving  inward.  After  that  time  or  r  >  0.24,  the  peak  is 

stable,  indicating  energy  balance  at  that  point,  that  is,  heat  supplied  by 
the  laser  pulse  merely  compensates  that  dissipated  by  diffusion.  The 
interface  receding  speed  increases  with  time  first,  reaches  its  maximum  and 
then  decreases . 
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If  the  triangular  pulse  shape  is  such  that  the  intensity  drops  off 
quickly  after  its  peak,  the  subsurface  temperature  peak  disappears.  As  can 
be  seen  in  Fig.  3(b)  and  Table  1(d),  the  temperature  decreases  monotonically 
with  the  depth  into  the  material,  and  consequently  the  internal  explosion 
can  not  occur. 

Since  the  shape  of  laser  pulses  in  practical  uses  is  more  Gaussian- 
like,  we  now  consider  a  pulse  with  intensity 

I(t)  “  Imaxexp(_t2/2a2)'  (8) 

where  a  is  a  parameter  characterizing  the  width  of  the  pulse,  and  1^^  is 

the  peak  intensity.  Numerical  results  for  a  YBCO  target  with  a  pulse  of  -30 
ns  long  are  shown  in  Fig.  4.  The  curves  in  Fig.  4(a)  illustrate  clearly  the 
subsurface  overheating  at  different  instants  as  long  as  the  pulse  intensity 
is  sufficiently  strong.  This  implies  that  internal  explosion  occurs  in  the 
laser  ablation  process  if  the  pulse  intensity  is  strong  enough,  and  hence 
provides  a  theoretical  basis  for  a  hypothetic  explanation  of  the  observed 
19 

solid  particles.  The  formation  of  such  subsurface  overheating  spots  can 
be  controlled  by  reducing  the  laser  intensity. 

According  to  parameters  reported  in  Ref.  5,  the  maximum  laser  intensity 

8  2 

*max  esC*-maCec*  to  be  -0.7  x  10  W/cm  ,  which  falls  within  the  range  of 

intensity  for  the  appearance  of  subsurface  overheating  in  Fig.  4(a).  The 
temperature  peaks  decrease  as  the  pulse  intensity  decreases.  They 

8  2 

eventually  disappear  for  Imax  <  0.13  x  10  W/cm  ,  which  may  be  regarded  as  a 

threshold  intensity  for  the  subsurface  overheating  under  experimental 
conditions  described  by  the  same  parameters. 
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Corresponding  to  every  instant  in  Fig.  4,  the  target  surface  receding 
velocities  for  all  the  cases  can  be  calculated  and  are  listed  in  Table  1. 

By  comparison,  we  find  that  the  receding  velocity  and  the  pulse  intensity 
reach  their  peaks  almost  simultaneously  for  high  laser  intensities.  On  the 
other  hand,  the  laser  intensity  reaches  its  peak  before  the  receding 
velocity  for  weaker  pulses,  as  expected. 

In  order  to  avoid  the  undesirable  products  of  internal  explosion,  the 
pulse  intensity  has  to  be  reduced.  Excessive  reduction  of  the  laser 
intensity,  however,  is  not  desirable  either,  because  the  evaporation  rate 
may  become  too  low  to  be  practical.  A  good  choice  is  probably  just  above 
the  threshold  intensity,  as  found  in  Fig.  4(c)  for  the  present  case.  It 
should  be  emphasized  from  the  above  discussions  that  the  pulse  shape  can 
also  make  an  important  difference,  especially  when  the  intensity  is  strong. 
In  Fig.  5,  we  compare  a  sine  wave  with  the  Gaussian  of  the  same  peak 
intensity.  The  difference  is  not  very  significant  because  Imav  is  not 

strong.  In  conclusion,  we  find  that  the  most  efficient  deposition  requires 
pulse  shapes  with  a  quickly  diminishing  tail  after  its  peak,  such  as  a 
narrow  Gaussian  pulse. 
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Figure  captions 

1.  Laser  pulse  shapes  considered  in  this  paper.  The  intensity  and 

8  2 

duration  are  measured  in  the  numerical  work  by  1.27  x  10  W/cm  and  30  ns, 
respectively,  (a)  Rectangular  pulse,  (b)  Triangular  pulse,  (c)  Gaussian 
pulse,  (d)  Sinuosidal  pulse. 

2.  Temperature  profile  as  a  function  of  the  depth  x-X  (in  60  A)  into  the 
kat  material  at  different  instants  (t  is  in  the  unit  of  30  ns).  The  switch- 
on  time  of  the  laser  pulse  is  taken  to  be  t  -  0.  Curve  1  is  for  the  moment 
when  the  surface  temperature  just  reaches  Tv;  2  is  for  t  -  0.08;  3  is  for  t 

-  0.16;  4  is  for  t  -  0.24;  5  is  for  t  -  0.32;  and  6  is  for  t  -  0.40.  (a) 
Triangular  pulse  with  a  -  1,  c  -  0.4  and  d  -  0.  (b)  Rectangular  pulse  with 

a  -  1  and  b  -  0.4. 

3.  Subsurface  temperature  versus  the  depth  x-X  at  different  instants  under 
irradiation  of  a  triangular  pulse  with  a  -  0.28,  c  -  0.028  and  (a)  d  -  1, 

(b)  d  -  0.17.  Otherwise,  all  curves  are  the  same  in  Fig.  2. 

4.  Subsurface  temperature  versus  the  depth  x-X  at  different  instants  for  a 
Gaussian  pulse  with  a  -  0.25  and  b  -  1.  Curve  1  is  for  the  moment  when  the 
surface  just  reaches  Tv;  2  is  for  t  -  0.2;  3  is  for  t  -  0.35;  4  is  for  t  - 

0.5;  5  is  for  t  -  0.65;  and  6  is  for  t  -  0.8. 

5.  Subsurface  temperature  versus  the  depth  x-X  at  different  instants  for  a 
(a)  Gaussian  and  (b)  sinusoidal  pulse,  both  with  a  -  0.2  and  b  -  1.  Curve  1 
is  for  the  moment  when  the  surface  temperature  reaches  Tv;  2  is  for  t  - 

0.35;  3  is  for  t  -  0.5;  and  4  is  for  t  -  0.65. 
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Table  1:  Interface  receding  velocity  at  different  instants  for  laser  pulse 


shapes  given  in  Fig.  1. 
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0.40 

1.32 

1.39 

19.7 

0 

16 


(b) 


at 


